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Abstract 

We calculate the pressure, entropy density, trace anomaly and speed of sound of the gluon plasma using the dilaton 
potential of Ref. [ 1] in the dilaton-gravity theory of AdS/QCD. The finite temperature observables are calculated from 
the Black Hole solutions of the Einstein equations, and using the Bekenstein-Hawking equality of the entropy with the 
area of the horizon. Renormalization is well defined, because the T = theory has asymptotic freedom. Comparison 
with lattice simulations is made. 
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1. Introduction 

The duality of string theory with 5 -dimensional gravity 
is nowadays a powerful tool to study the strong coupling 
properties of gauges theories and in particular of QCD, 
either at zero or finite temperature. In conformal AdS$ 
the metric is well known. It has a horizon in the bulk 
space at rj = ^- where {3 — 1/T and £ is the size of 
the AdS-space. Conformal solutions for entropy scale 
like s cc T 3 , since the 3-dimensional area of the hori- 
zon is given as A cc lip-pj. Promising solutions of this 
conformal theory have been proposed to the problem 
of viscosity rj J2[] with a small constant value for rj/s. 
Top to bottom approaches based on the conformal SYM 
with fermions have been investigating the chiral phase 
transition and problems at finite density yl |4J] . 

To extend this duality to SU(N C ) Yang-Mills theory, 
one of the first tasks is to control the breaking of con- 
formal invariance. In this work we will studied a 5D 
dilaton-gravity model introduced in [5], with the action 



S = 



— f 



(fx S-HK . 



The last term is the Gibbons-Hawking term, with K be- 
ing the extrinsic curvature and H the induced metric on 
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the boundary dM. The dilaton potential V(<p) is related 
to the /3-function in a one-to-one relation. While the 
yS-function is well known at high energies, there is no 
general consensus about its IR behavior. We assume the 
following parameterization 
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bio. + 



-#> or 
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-a I a 



(2) 



which was proposed in Ref. [1] for the computation 
of the heavy QQ potential at zero temperature within 
this dilaton-gravity model. This parameterization has 
the standard behavior in the UV region limit /3(a) 
-/3oa 2 - /3\a 3 + ... for a <K a, and generates confine- 
ment in the IR region, where it behaves as /3(a) -b^a 
for a » a. The optimum values to reproduce the QQ 
potential at T = are HI] 

b 2 = 2.3 , a = 0.45 , I = 4.389G6V- 1 . (3) 

This parametrization also allows to obtain the running 
QCD-coupling a(E) in the MS-scheme from the run- 
ning of the dilaton <p(z) in the bulk by mapping the bulk 
coordinate z to the energy scale in Ref. [1]. Small val- 
ues of the bulk coordinate z correspond to large energies 
and large values of z describe the infrared physics. 

2. Equations of motion 

The dilaton-gravity model of Eq. ([D has two different 
types of solutions for the metric. The thermal gas so- 
lution corresponds to the confined phase, and the black 
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hole solution characterizes the deconfined phase and it 
has a horizon localized in the bulk coordinate at z = Zh 
similar to the situation in 4-dim gravity [3]. We have 
studied in Ref. |6] the heavy QQ free energy in the con- 
fined phase of QCD by using a thermal gas metric. The 
black hole metric takes the form 

dz 2 



ds z 



b\z)\f(z)dT 2 



+ dxidx + 



(4) 



where z = l 2 jr and f(z) has the properties 

f(z h ) = 0, f'{ Zh ) = -4nT. (5) 

The second expression in Eq. (O follows from regu- 
larity of the metric at the horizon. This solution only 
exists at high enough temperatures [5]. The classical 
equations of motion corresponding to the BH solutoin 
read SI2I 
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where the superpotential is defined as W(z) = 
-9b'(z)/(4b 2 (z)). 

3. Thermodynamics 

The phase transition from the confined glueball gas 
to the deconfined gluon plasma can be qualitatively un- 
derstood as follows. At low temperatures the infrared 
physics is determined by the thermal gas metric with 
/ = 1. At a minimal temperature the black hole so- 
lution appears with a horizon at Zh and / + 1 and 
becomes the preferred solution at the phase transition 
which is first order. The phase transition temperature 
is fully fixed by the scales of the theory at zero tem- 
perature. The position of the black hole horizon inter- 
venes between the confinement region in bulk space and 
determines from then on the physics. Let us describe 
our new results starting with the highest temperatures, 
where expansions in the thermal coupling 07, = a(zh) 
are possible. Starting from the ultraviolet expansion of 
the /^-function, one can solve analytically the equations 
of motion in the UV. In doing that, one arrives at the 
following expansion for the metric factor up to 
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1 Analytical expressions up to 0(a 3 ) will be presented in Ref. 
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Figure 1: Pressure over T 4 as a function of T (in units of T c ). We 
show as red points the recent lattice data for SU(3) given in Ref. |9] 
for N% X N T = 20 3 X 5 , and as blue points the lattice data from Ref. flOfl 
for iV 3 X N T = 32 3 X 8. The colored curves represent the analytical 
result from the holographic model as an expansion in powers of a;,, 
c.f. Eq. 1141 . The order 0(a^) is computed in Ref. | 8]. The black 
solid line refers to the full numerical result. We have used the value 
of G5 quoted in Eq. )16t . 
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Figure 2: Entropy density over T 3 as a function of T (in units of T c ). 
See Fig.f^for details. 



The temperature dependence of the horizon is computed 
from the second expression in Eq. ©. It reads 



1 



l + ^a 2 h+ 0(al) 



(11) 



The entropy density follows from the Bekenstein- 
Hawking entropy formula which establishes the propor- 
tionality between the entropy and the area of the event 
horizon of the black hole. Up to Oia 2 ^ it reads 
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model reaches the Stefan-Boltzmann limit much slower 
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Figure 3: Trace anomaly density (e - 3p)/T 4 as a function of T (in 
units of 7V). See Fig.fT]for details. 



To reach this expression one has to evaluate Eq. cTTOt at 
the horizon, with Zh given by Eq. ( fTTT ). Using 



dT 



-Pwl-fral+0{ a A h ), 



(13) 



it is easy to derive the pressure from Eq. ( TTZl l. as it is 
related to the entropy density by s — dpjdT . One gets 



J3/3 



P(T) _ 

T 4 16G 5 



1-^0^ + 1(4^-3^)^ 



.(14) 



From Eqs. ( fl2l i and (fl4b one can immediately compute 
the energy density e - Ts - p, and the trace anomaly 
(e-3p)/T 4 = s/T 3 -4p/T 4 . 

In order to obtain the thermodynamic functions for 
all temperatures, numerical solutions of the equations 
of motion have been obtained §]. We show in Figs.Q] 
|2]and[3]the pressure, entropy density and trace anamoly 
for pure gluodynamics with N c - 3. The lattice data are 
taken from Refs. Hd. We also plot the analytical re- 
sults of Eqs. (fl2l i and (fl4l . It is noteworthy and visible 
in the figures that the expansion in terms of a/, converge 
quite rapidly. Even at temperatures T - l.5T c , two or- 
ders of a-/, give a sufficient approximation, quite oppo- 
site to the a(T) expansion in conventional high temper- 
ature QCD perturbation theory (pQCD). 

The question arises: which value of the gravitational 
constant G5 should we choose in the 5-dim gravity ac- 
tion? In principle, G5 can be chosen to reproduce the 
Stefan-Boltzmann limit at high temperatures, i.e. 



i6Gr 
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If one believes that the high temperature region T ~ 
(100r t . - 1000r c ) is a perturbative regime which can 
be described by the perturbative /J-function, then this 



than what lattice data suggest [11]. This is easy to see, 
because the coefficient of (9(a/,) in the expansion of the 
pressure in the holographic model (p^' s = P), Eq. (TPfl i. 
is a factor two larger than the corresponding one in 
pQCD (p\ QCD = £) [12]. As a consequence, the value 
of G™ given by Eq. (fT~5b leads to lower values of the 
thermodynamic quantities for all temperatures in com- 
parison with lattice data, also in the regime close to T c . 
To reproduce lattice data in the regime T c < T < 5T C 
one must use a value of G5 which is a factor 1 .46 smaller 
than G~, i.e. 



— = 1.46 x — 

G 5 G? 



(16) 



spoiling the Stefan-Boltzmann limit at high tempera- 
tures. Note, however, that lattice data for pressure and 
entropy density taken from Refs. ||9J] and II 1 OTI are not 
consistent each other at high temperatures, and it is 
hard to believe that both computations fulfill the Stefan- 
Boltzmann's law. This discrepancy introduces an error 
in G5 of the order of 10%, which in either case is not 
enough to explain the factor in Eq. (fTBI l. The gravity 
model seems to be hardly consistent to reproduce at the 
same time lattice data at very high temperatures, and 
close to the phase transition. In this sense, there exists 
the possibility that the ideal gas limit doesn't correspond 
to the limit of the black hole gravity theory at high tem- 
peratures. A natural question arises: is it possible that 
the gravity theory allows more degrees of freedom at 
high temperatures? We will further address this prob- 
lem, and analyze possible solutions |8[]. 

We plot in the figures as a black continuus line the nu- 
merical computation of thermodynamic quantities using 
for G5 the value quoted in Eq. ( TT6b . The trace anomaly 
in Fig.[3]shows the characteristic decrease towards zero 
for high temperatures expected from pQCD. Neat T c , 
for p, e and s in general all the AdS/QCD curves in- 
crease too slowly. 

The latent heat is defined as the energy density at T c , 
Li, = e(T c ). Another way to choose G5 would be to 
reproduce the value of L/, given by lattice simulations. 
In fact, with the value of G5 quoted in Eq. ( TTBT l one gets 
Lh/(T 4 ■ (N* - 1)) = 0.309, which is in good agreement 
with the result from lattice, 0.33(3) J9J]. 

Finally we can study the speed of sound c s which is 
independent of the normalization factor G5. From the 
specific heat per unit volume c v = and the entropy 
density s, the speed of sound writes 
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Figure 4: Speed of sound squared as a function of T (in units of T c ). 
We show as points the lattice data for SU(3) taken from Ref. [10] for 
N^. X N T = 32 3 X 8. Colored curves represent the analytical result of 
Eq. )17t and the black solid line refers to the full numerical computa- 
tion. The order 0(afy is presented in Ref. 1 8]. 

In the r.h.s. we write the result of the computation of this 
quantity in the ultraviolet from Eqs. ( fl2b and ( fl4l >. We 
show in Fig. |4] the speed of sound computed with the 
holographic model, and compared with the lattice data 
of Ref. Midi . We also show the analytical ultraviolet ap- 
proximation given by Eq. ( fTTb . Since c] is close to 1 /3 
in the calculation, we see that we have massless excita- 
tions in the plasma in the range 2T C < T < 5T C . 

4. Conclusions 

We study in this paper the thermodynamics of the five 
dimensional dilaton-gravity model by using a parame- 
terization for the dilaton potential which was quite suc- 
cessful to reproduce the heavy QQ potential at zero tem- 
perature HI. We compute analytical expressions for the 
pressure, entropy density and speed of sound, as an ex- 
pansion in powers of the running coupling. This ex- 
pansion turns out to converge quite rapidly even at tem- 
peratures T ^ l.5T c , quite opposite to the conventional 
QCD perturbation theory at high temperature. The grav- 
ity model with the dilaton potential of Ut] cannot repro- 
duce both the low temperature T x T c and the Stefan- 
Boltzmann limit at very high temperatures. This finding 
differs from the result obtained with a different dilaton 
potential shown in Ref. II 1 3fl . Both potentials are based 
on /^-functions which agree in leading order. The un- 
derlying /3-function of Ref. 11311 has a coefficient /fee? 4 
which is rather large. Consequently the physics based 
on this model shows strong nonperturbative features al- 
ready at very small coupling. It is well known that the 
coefficient @2 is scheme dependent and it is possible to 
set the physical scale for the thermodynamics also in 



this model from the QQ-potential. The mapping, how- 
ever, between the a obtained from the gravity theory 
and has not been achieved. The thermal couplings 
a(zh) in the gravity model of Ref. [13] are very much 
smaller than the a(T) in the corresponding QCD-lattice 
simulations, thereby making the calculation compatible 
with the asymptotic Stefan-Boltzmann pressure. The 
swift change from a perturbative to a nonperturbative (3- 
function facilitates the steep rise of the thermodynamic 
functions at low temperatures. 

The description of Ref. [1] does well for all quan- 
tities which are calculated in the string framework, like 
the free energy of a QQ- pair. In a forthcoming work J^l] 
we will demonstrate the agreement between the compu- 
tation of the free energy from the Bekenstein-Hawking 
entropy formula presented here, and the method fol- 
lowed in Ref. |5[] based on the regularization of the 
Einstein-Hilbert action of Eq. ([1). We will also extend 
the computation to other thermodynamic observables, 
like the Polyakov loop and the spatial string tension. 
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